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C^ . Abstract 



The opportunity to observe color transparency (CT) is determined by how 
rapidly a small-sized hadronic wave packet expands. Here we use SU(2) lattice 
gauge theory with Wilson fermions in the quenched approximation to investigate 
the expansion. The wave packet is modeled by a point hadronic source, often used 
as an interpolating field in lattice calculations. The procedure is to determine the 
Euclidean time (t), pion channel, Bethe-Salpeter amplitude ^(r, t), and then evalu- 
ate 6^(t) = / d^r^(r, t)r^sm^0^7r(^)- This quantity represents the soft interaction 
of a small-sized wave packet with a pion. The time dependence of 6^(t) is fit as a 
superposition of three states, which is found sufficient to reproduce a reduced size 
wave packet. Using this superposition allows us to make the analytic continuation 
required to study the wave packet expansion in real time. We find that the matrix 
elements of the soft interaction b^ between the excited and ground state decrease 
rapidly with the energy of the excited state. 
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1 INTRODUCTION 

Color transparency (CT) can be defined as the absence of or reduced final state interac- 
tion of a nucleon witli tlie nucleus in semiexclusive liigli momentum transfer processes. 
Examples are the (e,e'p) and (p,pp) reactions occuring on nuclear targets, in which the 
detected protons have enough energy to ensure that no pions are produced. For color 
transparency to be observed three conditions have to be met |jl|. 

(i) A small wave packet is formed in a high momentum transfer reaction. This wave 
packet is sometimes dubbed a point like configurations(PLC). 

(ii) A small wave packet interacts weakly with the nucleus. 

(iii) The wave packet escapes the nucleus before expanding. 

These conditions are true in the perturbative regime p|, 0, but have to be tested in 
the nonperturbative regime. Lattice QCD is a natural tool to use in this investigation. 
It can be used to investigate the form of the hadron-hadron interaction [Q, to probe the 
wave packet (hopefully PLC) produced in a high momentum transfer reaction, and to 
examine the expansion of the PLC. 

The separation of the interaction into the hard high momentum transfer part and the 
soft final state interaction part is crucial for CT. To be specific consider (e,e'p) reaction. 
In this case the hard process is the absorption of the virtual photon by a nucleon in the 
nucleus. An additional hard interaction is expected to be suppressed by the experimental 
kinematics, and the final state interaction of the PLC with the nucleus is of a soft low- 
momentum-transfer nature. There are two amplitudes to be compared. The first one 
is for the proton to escape the nucleus without interaction. The second one is for the 
proton to be detected after the PLC is scattered by the soft interaction with the nuclear 
medium. CT can be obtained if the ratio of the second amplitude to the first one vanishes 
in the limit oi Q —^ oo. 

Ideally one would like to reconstruct on the lattice the process described above: to 
create a wave packet by acting with the electromagnetic current operator on the proton, 
and to evaluate the amplitude of this wave packet to be transformed into a proton by 



some realistic nucleon-PLC interaction. In this paper we will not take on this challenging 
and complex problem. We concentrate on one aspect of the question - the expansion of 
a model small size wave packet. To further simplify our task we consider two-quark pion 
instead of three-quark proton. 

We start by defining the necessary lattice constructions in Sec. 2. A point source, one 
of the conventionally used interpolating fields, is used to generate the small size wave 
packet - PLC. A pion is detected in the final state after the PLC experiences the soft 
interaction. A model form of the soft interaction b^ [Q is used, with b representing the 
transverse separation of the quarks. The wave packet undergoes Euclidean time evolution. 
The quantity &^(t), which measures the strength of the soft interaction of the evolving 
PLC with the nuclear medium, is calculated as a function of the Euclidean time t. The 
analytical continuation to the Minkowsky time r is then performed. In Sec. 3 the lattice 
details are given. The purpose of using lattice QCD is to obtain matrix elements from the 
first principle calculations. Unfortunately it entails the Euclidean time evolution. Our 
procedure is to treat the wave packet as a coherent sum of physical states, so the real 
problem is how to detect as many states as possible before they decay away. These states 
are extracted by a many pole fit of the lattice results. This procedure is complicated and 
controversial. We discuss it for the two-point correlation function in Sec. 4. To perform 
the analytical continuation the b'^{t) is fit with the sum of three exponentials. The details 
of the b'^{t) three pole fit are given in Sec. 5. The results of the analytical continuation and 
their discussion are presented in Sec. 6. The relative merits of using a full or a diagonal 
covariance matrix to fit observables is discussed in the Appendix. 

2 FORMALISM 

We want to consider the expansion of a small size wave packet. To form such a wave 
packet we take advantage of the point interpolating field 0(x, t), which has the following 



form in the pseudoscalar channel 

0(x, t) = d{x, t)75M(x, t) (1) 

By acting on the vacuum 0(x, t) creates an object similar to a PLC, since the quark 
and the antiquark are created at the same site. This wave packet can be expanded in a 
complete set of zero momentum states |z/ >: 

^0(x,O)|O> => Y.(j dp)c,\u>=\PLC> (2) 

where the summation over x projects out the zero momentum states. The coefficient Cy^ 



c, =<i/|0(O,O)|O>/^2E,, (3) 

is the strength with which the state |i/ > is excited from the vacuum, \p2El is a normal- 
ization factor. We have to note that even though the interpolating field creates a quark 
and an antiquark at the same site, there is a non-zero probability to detect them sepa- 
rated by some distance at the same time they are created. This occurs because fermion 
propagators do not vanish outside of the light cone. 

If we transport the \PLC > with the Euclidean transfer matrix exp{—Ht) and then 
contract it with itself we obtain the two-point correlation function G2{t): 

X 

= Y.i[du)\cy\'e~^^' (4) 

6*2 (^) plays a complementary role in our investigation, as will be seen below. 

Our main objective is to evaluate the strength of the soft interaction of the expanding 
PLC with the nucleus. This interaction can be described (0 and references therein) by 
the operator b"^ with b representing the transverse separation of the quarks in the PLC Q 

b\t) = -< TTleH^e-'^'lPLC > . (5) 

Cl 



^For the real PLC the longitudinal direction is given by the momentum transfer in the hard inter- 
action. Unlike real PLC our wave packet is spherically symmetrical so we arbitrarily choose the z— 
direction to be longitudinal. 



Here we normalize it by ci =< 7r|PLC >, which represents the amphtude of the detecting 
the pion without the soft interaction. To find the matrix element ( |^) we need to know 
the wave functions of the PLC and the pion. There is no unique way to define the 
wave functions on the lattice. We use the gauge invariant formulation of equal time 
Bethe-Salpeter amplitude \l/(r, t) 

^(r, t) = ^ < 0|S(x, t)75f/(x ^ X + r, t)u(x + r, t)0(O, 0) |0 > . (6) 

X 

The Bethe-Salpeter amplitude in the gauge invariant formulation is known [0 to 
underestimate the spatial extent of the pion wave function. But it does not suffer from 
finite lattice size errors and is hoped to give a qualitatively correct estimate of the effects 
under consideration. 

In a fashion similar to G2{t) ^^(r,^) can be expanded in a complete set of states 

^{r,t)=Y.{fdu)c,M^)e-^'\ (7) 

where 



V-.W =< 0|rf(0,0)75?7(r,t)u(r,t)|z/ > / pE, (8) 

is the Bethe-Salpeter amplitude of the state |z/ >. 

\l/(r, t) is measured on the lattice and lP'{t) is calculated according to 

b'^it) = f dh^{r,t)ilj^{r,t){rcosey. (9) 

The ground state wave function ■?/'7r(r,t) is obtained from \I^(r, t) in the limit t -^ oo. 
6^(t) has the following spectral representation 

h\t)=Y.{ldvf-^hle-^^^-^^^\ (10) 

where 

bl = J d^rtP,{r,t)4j^{r){rcosef. (11) 

The spectral representation ( p!0| ) can be used to perform the transition from the 
Euclidean time to the Minkowsky time. If we knew the E^^s, c^^s, and b^^s for all the 



contributing states, then a mere substitution t ^ ir would give us &^(t), which describes 
the expansion of the PLC in the real time. This is certainly too much to ask for. Since we 
work with the Euclidean time most of the states decay away too quickly to be detected. 
We assume that the spectral density can be represented as a sum of several sharp 
poles. With this assumption the goal is to try to determine the E^^s, Cj^'s, and b^'s for 
as many states as possible. We perform the many pole fit of &^(t) with the fit function 
Mt): 

Np 

/,(t) = ^a.e-^"*, (12) 

n 

where Np is the number of the poles. However, there are several sticky points here. First, 



there has been some concern expressed in the literature [^ about the validity of such 
a many pole fit. Second, it is not a priori obvious that those few states we are able to 
extract will be enough to form a small or at least reduced size wave packet. In Sec. ^ 
and Sec. ^ we show that a simultaneous three pole fit of b'^{t) and 6*2 (^) can be reliably 
performed to yield the parameters for the three lowest states. These three states do form 
a reduced size wave packet, whose expansion is considered in Sec. ^ But first we give a 
summary of the lattice calculations in the next section. 

3 LATTICE DETAILS 

The calculations are performed for SU{2) gauge theory with the coupling constant (3 = 
2.5. SU{2) is chosen to increase the efficiency of the calculation and improve statistics. 
The size of the lattice is 12^ x 24. 60 quenched gauge field configurations are generated 
using Metropolis method with overrelaxation 0. The first configuration is selected after 
2000 thermalization sweeps and all the consecutive ones after 1000 sweeps. 

The Wilson propagators are calculated for three values of the hopping parameter k, = 
0.146,0.148,0.149. Periodic boundary condition in the spatial directions and Dirichlet 
boundary conditions in the time direction are used. A standard procedure is used 
to construct the two-point correlation function and the Bathe-Salpeter amplitude out of 
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the quark propagators. 

The lattice spacing is determined from the p mass [§ to be a = 0.09 ± 0.012 fm, with 
is in quantitive agreement with the results 0] obtained from the string tension. 

The quark mass mq is conventionally associated with the hopping parameter k: rUg = 
(1/2/t: — l/2/tcr)a^^. Extrapolation to the limit rrij^ = yields K^r = 0.151(1). With this 
Kcr we obtain three values of the quark mass rriq = 275, 174, 124 Mev. 

The covariance matrix fit is performed to determine the parameters needed. A de- 
tailed description of this procedure is given in Sec. 0. All statistical errors are estimated 



by the single elimination jacknife |10|. The point source is placed at the time slice t = 5 



and the fitting is performed over the range t = Q through t = 20. 



4 MANY POLE FIT 

A many pole fit involves many problems, some of which are mentioned in [|l^, |lT|. To ad- 



dress these problems we perform a detailed analysis of the fit of the two-point correlation 
function G2(t). The fit function fcit) 

Np 

/gW = E Cne-^"* (13) 

n=l 

corresponds to A^p-pole fit. 

We would like to emphasize that our goal here is to extract as many states as possible 
to be able to form a small size wave packet. This is contrary to the intention of the 
majority of the papers in the field, where the isolation of one state, usually the ground 
state, is pursued. 

A standard approach in determining parameters c^'s and E^s is to minimize the 
goodness-of-fit statistic x^ ■ 

X' = E(G'2(t) - fGit))C{t,tr\G2it') - fait') (14) 



The jacknife method is used to calculate the covariance matrix C{t,t') 



N, 



Cit,t') = ^^^3Y E(G^2 W - G,{mG['\t') - G,it% (15) 



where Nc is the number of configurations, G2 (t) are found for a subset of A''^ — 1 config- 
urations with k'th configurations omitted. 

Simulated annealing is our choice of the minimization technique. The full covariance 
matrix is used in the minimization(correlated x^ fi^)- There have been studies [|l^, |17 



to show that for a small data sample use of the full covariance matrix is unreliable. 



In |TT|] a use of several well defined eigenmodes of the covariant matrix is proposed. 
The use of a full covariant matrix versus the diagonal part thereof poses an interesting 
problem. We address this problem in the Appendix, where we consider a simple example 
amenable to the analytical treatment. The conclusion is reached that a correlated fit 
gives better results if the correlations are appreciable. This is generally the case for lattice 
calculations, and it is true for our calculations. Moreover, we consider 60 configurations 
to be sufficient for the maximum of 9 degrees of freedom we have while fitting 6*2 (t). 

The fitting is performed over a time range extending from t first to tiast- We choose 
tiast = 20 to exclude the boundary effects and t first is varied from 6 to 14. The number 
of states used in a fit is determined by the t first- One of the problems with a many pole 
fit [jl2| is that the result depends on the t first used for the fit, and it is not a priori clear 



how many states should be included in a fit. 

The strategy developed to tackle this problem is as follows. We start with t first for 
which only the ground state contribution is significant, and one-pole fit is performed. As 
t first decreases more states with higher masses come into play. To decide between two 
fits with different number of poles for a particular t first the preference is given to a fit 
with smaller value of x^ per degree of freedom. The results for the En's for the hopping 
parameter k = 0.149 are shown in Fig. 1. For tfirst = 12 -^ 14 one state is enough, for 
t first = 9-^11 two states had to be included, and for tfirst = 6^8 the G2{t) is best fit 
with three states. The results for the other two hopping parameters are similar. The 
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values of the x^ P^^ degree of freedom show that the fit for the longest time range is the 
best( Table 1). 

The pattern in the results for En is obvious. The statistical error is the greatest 
when a new state is included in the fit, it then decreases and jumps back, when another 
state is included. To explain this pattern the following plot is made. The individual 
contributions a„e~^"* of the states to the two-point correlation function G2{t) are shown 
in Fig. 2. They are compared with the statistical error on the 6*2 (t). Two time slices 
ifirst = 12 and tji^st = 9 are important, since at those times the second and third 
states, go above the statistical noise. Those are precisely the time slices for which the 
contribution of these states is included in the fit. 

A simple rule follows: a state should be included if its contribution is greater than the 
statistical error; it should not be included if its contribution is less than the statistical 
error. The best results of the fitting are obtained for the longest time span of the highest 
energy state used in the fitting. An interesting consequence is that, in general, the 
longest time range is not necessarily the best for determining the parameters. Consider 
the following example. Suppose that in Fig. 1 the results of the three pole fit started 
at tfi-fst = 8 (i.e. the source was at t = 7 and the results for tfirst = 6 and 7 were not 
present). In that case a reliable determination of the third state's parameters would be 
impossible; furthermore we would be forced to use the shorter time range with tfirst = 9 
to determine the parameters of the two lowest states. 

Another interesting observation concerns the fact that x^ has many minima in the 
multi-dimensional space of c„'s and E^s. The "real" minimum is not always the global 
one. By the "real" we mean the one we believe to be the best fit. In our case the best fit 
is obtained for the tfirst = 6. We ran the minimization for tfirst = 7 and tfirst = 8 and 
found local minima very close to the "real" one. The gain in x^ compared to the global 
minimum is only ^0.1 per degree of freedom. 

To support this argument the following simulation is performed. We work our way 



backwards. Four sets of c„ and En are chosen. A "correlation function" G{t) is calculated: 

G{t)= E 5nexp(-E„t). (16) 

n=l,4 

Then statistical noise is added and 60 "configurations" are generated. The subsequent 
fitting procedure yield results that have all the basic features of the fitting of the lattice 
data(Fig. 3). The input values of the parameters of three states are recovered (within 
the error bars). The fourth state is intentionally chosen to decay below the statistical 
noise by t = 6. In this simulation we know what the correct values of the parameters 
are, and we can verify that for tfirst = 7 and tfirst = 8 there are local minima very close 
to these correct values. 

We conclude that for the given choice of interpolating field and lattice spacing three 
pole fit of the two-point correlation function can be performed and the values of the 
excitation strength c„ and the mass En reliably determined. These values are presented 
in Table 2. 

5 b\t) MANY POLE FIT 

The determination of the c„'s and EnS is interesting, but it is does not constitute the 
main objective of this paper. What does is the strength of the soft interaction &^(t), 
which we fit with the Eq.( |1^). As was indicated in the introduction, a fit of &^(t) alone 
does not give a satisfactory result. To understand it, note that a coefficient a„ in Eq.( |T^ 
is proportional to the excitation strength c„ and the soft interaction transition matrix 
element 6^ of the state \n >. For the three states, that we are able to extract from the 
two-point correlation function G2(t), the excitation strength c„ grows with n. (Physically 
it means that their wave functions have a better overlap with the wave function of the 
initially formed wave packet.) Unlike the coefficients Cn h"^ turned out to decrease rapidly 
with n. As a consequence the contribution to the 6^(t) of the exited states disappears 
below the statistical noise much sooner than for the G2{t). This forces us to fit 6^(i) 
together with 6*2 (i). 
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We would like to note here that the procedure we use to compute 6^(t) is akin to the 



smearing technique used in [^, |IJ] to isolate the contribution of meson and baryon ground 
states. However, let us stress it again, this is not our intention. This is an "accident", 
but an "accident" with some physical consequences, which are discussed below. 

The simultaneous fit of &^(t) and 6*2 (^) is performed for the longest time range possible 
{t first = 6). The results are shown in the Table 3. 6^ goes down by a factor of 2 with 
every excited state. Another important feature is that 6| is negative. This allows for the 
formation of a reduced size wave packet. 

Let's discuss the physical meaning of these results. Do they have any relevance for 
color transparency? Again consider (e,e'p) process. In this case the coefficients c^, defining 
the PLC are the elastic and inelastic form factors of the ground state. In general we can 
expect these coefficients to exibit a resonant behavior - to have a peak for a certain excited 
state. But the soft interaction supresses the contribution of the higher excited states, as 
biy decreases rapidly with u [|l^. This suppression plays a crucial role in achieving a slow 
expansion rate of the PLC. Our result turned out to have these basic features. The c„ 
grows with n, whereas 6„ decreases rapidly. The former is only a result of our choice of 
the model wave packet and in this sense is arbitrary. The latter, however, reflects the 
properties of hadronic states in lattice QCD, and to the extent that we believe lattice 
QCD, this is an indication of the possibility of the suppression of the contribution of the 
higher excited states by the soft interaction. 

6 MINKOWSKY TIME EXPANSION 

Having obtained the parameters En,Cn, h\ for the three lowest states, we substitute t —^ ir 
in Eq.( |l^) and plot the real part of &^(r) (Fig. 4). 

Re{b\r)) = J2 -blcos{E^ - E,)t (17) 

The band shown in Fig. 4 corresponds to the 60 sets of 59 conflgurations used in the 
jacknife evaluation. Each of the curves shows an expansion. We do not compute a set of 
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curves simply using the errors quoted in Table 2 and Table 3. Such a procedure would 
incorrectly ignore the correlations between the parameters En, Cn, and 6^. 

A look at Fig. 4 shows us that the wave packet under investigation is not really 
small. According to Fig. 4, 6^(0) ~ 3/4 bl, where bl =< 7r|6^|7r > is the ground state 
contibution. (In the Fig. 4 6^ is represented by the solid horizontal line.) This is because 
we could recover only three out of the many states comprising the wave packet at the 
source t = 5. Furthermore, even at t = 5 the wave packet is not quite point like. As was 
mentioned above , there is a non-zero probability for the quark and the antiquark to be 
separated by some distance at the time they are created by the point interpolating field. 

These features make it difficult to define one number r^xp to typify the time for the 
wave packet to expand. One way is to define Texp as the time it takes for &^(t) to reach 
the ground state contribution bf. This expansion time is found to be ^ 0.03 fm, which 
is certainly much smaller than the root mean square radius of the pion ^ OA ira here 
0. However, this estimate is misleading, since the initial wave packet is not small. The 
initial wave packet size is determined by the coefficients c„ as well as b^- Here the set 
Cn determines the spectral density of the operator (f){x). These numbers are only related 
to color transparency physics if they also represent the elastic and inelastic form factors. 
Thus we are free to consider the consequences of using different values of Cn- Varying 
Cn we can obtain an upper limit for the expansion time for the given energies En and 
soft interaction matrix elements b^- This is achieved for C3 = and C2 < and equal to 

Texp ^~-' ^ im. 

The results for the other two hopping parameters are similar to the one shown in 
Fig. 4, but with bigger statistical error. For the purpose of our investigation we are not 
interested in taking the physical limit Ex —>■ ??i^. The point is that the expansion time 
is determined by the energy differences En — Ei as well as c„ and 6^. The coefficients Cn 
and bn are virtually independent of the quark mass. The lowest value of the energy gap 
^ This small number is the consequence of using the gauge invariant formulation of the Bethe-Salpeter 
wave function. 
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for a pion is much larger, than for a proton. By taking the physical limit to m.^ ^ we 
only increase the mass gap and decrease the expansion time. 

There is one more point to be made. The subject of our investigation is the expansion 
of a wave packet in its rest frame. If color transparency is to be observed, it will be in 
a high momentum transfer reaction, when the escaping wave packet is moving with a 
speed close to the speed of light. Time dilation will increase the expansion time. In fact, 
color transparency can be achieved for any small size wave packet, consisting of a finite 
number of states, as long as the momentum transfer is much greater than the energy 
of the states. A real PLC produced in a high momentum transfer reaction consists of 
an infinite number of states. How the energy of the important states changes with the 
momentum transfer will decide the fate of color transparency. 

7 CONCLUSION 

The expansion of a reduced size wave packet on the lattice is investigated. The wave 
packet is created with the pseudoscalar point interpolating field. The quantity represent- 
ing the strength of the soft interaction of the expanding wave packet with the nucleus is 
measured as a function of the Euclidean time. We treat the expanding wave packet as 
a coherent sum of physical states. Most of the states decay away before we can detect 
them. The three lowest energy states we are able to recover allow for the formation of 
a reduced size wave packet. With the parameters of these states we are able to do the 
transition to the Minkowsky time. A Minkowsky time expansion picture is obtained. 

What is the importance of our work and where do we go from here? For the first 
time three states - a ground state and two excited states- are obtained from the lattice 
calculations. (One may attempt to recover more excited states by using a smaller lattice 
spacing a.) The matrix elements of the soft interaction between the ground state and 
the two excited states are important for color transparency. These are found to decay 
rapidly with energy. This feature is an important component of color transparency. (It 
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would be useful to verify that similar results are obtained for baryons.) Furthermore, the 
present results show that a lattice evaluation of wave packet expansion is possible. 

In the present results, the initial wave packet is not small, and its relation to color 
transparency physics is unclear. Our next step will be to attack this problem by studying 
the wave packet (hopefully a PLC) created by the electromagnetic current operator. 

It is our pleasure to thank Steven Sharpe for useful discussions and encouragement. 
This work was supported in part by the U.S. Department of Energy. The lattice calcu- 
lations were performed on the UW Nuclear Theory DECstation 3000-600 AXP. 

A FULL VS. DIAGONAL COVARIANCE MATRIX 
FIT 

There have been several papers, in which the attention has been brought to the fact that 
the use of the full covariance matrix for the correlated data is unreliable for small data 
samples [ff6l IT^ . We use a simple example to show analytically that there are instances. 



when the use of the diagonal part of a covariance matrix to fit data (uncorrelated x 
fit) with a theoretical expression gives a better result than the full covariance matrix 
(correlated x^ fit), even in the limit of infinite number of samples. Note that when we 
use the diagonal part, we fit correlated data and simply ignore the correlations among 
them. 

Consider fitting results of the N sample measurements of A^ quantities x by a constant 



-^^ sample 



I. We assume that all Xj's 

^=Tr^ E ^? (18) 

-'■^sample a=l 

do reach the same constant value J in the limit Ngampie -^ oo and that the correct infinite 
Nsampie form of the covariance matrix Cik is known. 
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Then we can minimize x^; 

X^ = J2(xl-I)ak(xj:-I), (19) 

ik 

with respect to / and find I in two cases: I/uii, when the full covariance matrix C is 
used, and Idiag, when only the diagonal elements of the covariance matrix C are used. 

J-full — ,^ ^_i ; J-diag — ^^ -, ,^ ■ l,^UJ 

Both //„// and Idiag are unbiased estimators of the true value J, since in the limit 
^sample ~^ oo they both are equal to J . The preference should be given to the estimator 
with the smaller standard deviation. Results for /j„;; and Idiag will vary from one set of 
Nsampie Samples to another with the variance Cj = P — (J)^ equal to 

I 

Cifuu = C^c^k y Ci^^^^ = C^-Tj-y ■ (21) 

ik i ** 

From this general form it is clear that we can not expect C/ < Cj^^^ to be always true. 
Consider a situation, when the diagonal elements of the matrix C are much greater than 
nondiagonal: 

^«1; %^«l; ^^k, (22) 

^ii ^ii 

and we can expand C^^ to get for Cl^^^^ 



Ci,.u = {Y.^-Y.^c.>^r' (23) 

i ^tt i^k ^kk 



The uncorrelated fit yields better results in the case of prevailing positive correlations, 
the correlated fit is better in the case of prevailing negative correlations. But positive 
correlations do not necessarily render a correlated fit worse than uncorrelated one. Con- 
sider a simple example of A^ = 2, when the expressions for C/ and Cj^.^ become 
transparent: 



'-^ll'-^22 ~ '-^12 ^ _ '-^ll'-^2 

(^ I (~< oT^f Idiag ~ ~F< I r~< 



^ I full ^ \_ n on Idiag /^ \_ H \ J 
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Again, as in the case of small correlations, a correlated fit is better for negative correla- 
tions Ci2 < 0. But there is a region of big positive correlations 

C. > j|%^ (25) 

where a correlated fit gives better results. 

Another issue is using x^ to estimate how good a fit is. For the correlated fit the 
infinite N sample limit of Xfu// is A^ — 1, which is exactly 1 per degree of freedom. In the 
case of the uncorrelated fit x^ reaches this limit only if the data is genuinely uncorrelated 
{Cik = 0,i 7^ A;). When we ignore the correlations and use the diagonal part of the 
correlation matrix to obtain the result for J, they pop up in the expression for Xdiao- 

xL, = AT - 1 - E -7^/Y.iCu)-' (26) 

One therefore can not always use x^ to evaluate an uncorrelated fit. 
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Parameter ki = 0.146 K2 = 0.148 ^3 = 0.149 



Eia 



E2a 



E^a 



C\ 



C2 



C3 



0.49(1) 


0.39(4) 


0.33(5) 


1.2(1) 


1.1(2) 


1.0(2) 


2.2(1) 


2.2(2) 


2.1(2) 


0.39(3) 


0.36(6) 


0.35(6) 


1.1(2) 


1.1(2) 


1.0(3) 


2.3(1) 


2.4(1) 


2.4(6) 



Table 2: Energy E^ and excitation strength c„ of the three pole fit of the two-point 
correlation function 6*2 (t) 
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Figure 1: The energies of three states as a function of the first time shce t first used in 
the fitting, tiast = 20 for all the fits. The dash line shows the result for the longest time 
range. Hopping parameter k, = 0.149. 
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Figure 2: Individual contribution of the three states to the two-point correlation func- 
tion to be compared to the statistical error of the two-point correlation function 6*2 (t). 
Hopping parameter k, = 0.149. 
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Figure 3: The energies of three states for the simulated G{t). tiast = 20 for all the fits. 
The dash lines show the correct input values of the energies. 
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Figure 4: Expansion of the wave packet in real time. The sohd curve corresponds to 
the average over all 60 configurations. The dash curves are obtained for 60 sets of 59 
configurations that were used in the jacknife evaluation. The horisontal solid line is the 
ground state contribution bf. Hopping parameter j3 = 0.146. 
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Parameter ki = 0.146 K2 = 0.148 K3 = 0.149 

bj 0.08(2) 0.08(2) 0.08(2) 

bl 0.04(1) 0.05(2) 0.05(2) 

bl -0.02(1) -0.02(1) -0.015(13) 

Table 3: The soft interaction matrix elements 6„ =< nlfe^ll >. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9411030v2 
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http://arXiv.org/ps/hep-lat/9411030v2 
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http://arXiv.org/ps/hep-lat/9411030v2 
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